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^ : Abstract 
in 

' The preequilibrium (nucleon-in, nucleon-out) angular distributions of ^"^Al, 

^^Ni and ^^Zr have been analyzed in the energy region from 90 to 200 MeV in 
^ I terms of the Quantum Moleculear Dynamics (QMD) theory. First, we show 



that the present approach can reproduce the measured (p,xp') and (p,xn) 
^ I angular distributions leading to continuous final states without adj using any 

parameters. Second, we show the results of the detailed study of the preequi- 
librium reaction processes; the step-wise contribution to the angular distri- 
bution, comparison with the quantum-mechanical Feshbach-Kerman-Koonin 
theory, the effects of momentum distribution and surface refraction/reflection 
to the quasifree scattering. Finally, the present method was used to assess 
the importance of multiple preequilibrium particle emission as a function of 
projectile energy up to 1 GeV. 
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I. INTRODUCTION 



The nucleoli- induced nuclear reactions leading to continuum states at intermediate 
(Einc >100MeV) energy range are characterized by a reaction mechanism known as the pree- 
quilibrium process |jl|] . In this process, particle emissions take place from simple particle- hole 
configurations populated as a result of a sequence of nucleon-nucleon interactions before the 
statistical equilibrium is attained. The angular distribution of the particles emitted from 
this process has generally a smooth forward peaking which is intermediate in character be- 
tween the direct and compound nuclear processes. As the energy of the projectile increases, 
the number of particles emitted from the preequilibrium mechanism is increased and exceeds 
one (which is therefore called as the multiple preeqilibrium emission process, MPE process). 
At very high energy, the reaction is often referred to as the "spallation" reaction, in which 
the average multiplicity of ejectile exceeds several or larger. 

Study of the preeqilibrium nuclear reactions has been an active field since the pioneering 
work of Goldberger and Metropolis ||^ based on the cascade model, and of Griffin 
based on the exciton model. Various refinments on these approaches as well as new models 
both of semi-classical and quantum-mechanical followed (see for example Ref. 0). The semi- 
classcial models have been applied to analyze the energy spectra of preequilibrium particles 
on the outset. Later they have been improved to take account of the angular distributions 
of the preequilibrium process; the exciton model was improved to the generalized excition 
model and the geometry dependent hybrid model |§| . The cascade model has been able 
to calculate the angular distributions based on the Monte-Carlo technique. Furthermore, a 
semi-classical distorted wave theory was proposed by Luo and Kawai who have com- 

bined the concept of quantum distorted-wave and the cascade model. They have applied this 
theory to calculate the 1-step double-differential cross sections. Extension to the 2-step pro- 
cess was also formulated . Althouth these theories gave overall agreements with the data, 
there are still open questions which need further investigation for a better understanding 
of the preequilibrium reaction processes. For example, the backward angular distributions 
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calculated by the semi-classical theories are often considerably smaller than the measured 



values. Various conjectures have been made to account for this problem p,|T0|,[T2|-|T^] ; the 
refraction effects at the nuclear surface, quantum diffraction, high momentum component in 
the momentum distribution, multi-step effects, MPE, etc. So far, no simple answer seems to 
resolve this problem. The same problem persists even in the quantum-mechanical Feshbach- 
Kerman-Koonin (FKK) theory one and the only quantum- mechanical preequilibrium 
theory which is able to calculate the multi-step direct process up to any number of steps 
at present. In the FKK theory, furthermore, there are some other open problems, e.g., the 
transition between the unbound and bound states (P <^=^ Q transition) as studied recently 
by Watanabe et al. Jl^, and use of the normal and non- normal DWBA matrix elements in 



the calculation of multistep direct components [0 . On the contrary, the cascade model has 
a problem at both the very forward and backward angles, where the calculated values are 
noticeably smaller than the experimental data. Moreover, the number of particles emitted 
from the preequilibrium process is limited to only 1 or 2 in many preequilibrium theories 
proposed so far |T^|T^; an assumption which is questionable when the projectile energy 
becomes higher and higher. 

The purpose of this paper is to study the angular distributions and MPE process in 
the preequilibrium (nucleon-in, nucleon-out) reaction in terms of a reaction theory based on 
the molecular dynamics technique, the Quanum Moleculear Dynamics (QMD) p0|^2|. The 
QMD theory includes, in a self-consistent way, many important aspects in understanding 
the nucleon-induced reaction mechanisms at intermediate energy range, i.e., 1) the realistic 
momentum distribution of nucleons inside nuclei (including high-momentum component), 
2) entrance/exit channel refraction, 3) Coulomb deflection, 4) multistep process, 5) MPE 
process, 6) variation of the mean-field potential due to particle-hole excitation and par- 
ticle emission, 7) transition between unbound and bound states (P <^=^ Q transition 
in FKK theory), and 8) energy-dependent, anisotropic N-N elastic and inelastic scattering 
including the Pauli-bocking effect. These features make QMD a useful tool to study the 
nucleon-induced preequilibrium processes in a systematic manner as was first demonstrated 



by Peilert et al. |22]. We are willing to show calculations of better statistics to check its 



ability at very backward angles for energetic ejectiles which was not clear in Ref. ||22|. Fur- 
thermore, we will clarify the roles of the step-wise contributions, MPE contributions, the 
momentum distribution and surface refraction/reflection to understand the basic physics of 
the preequilibrium reactions. 

In this paper, we use the method as formulated in Ref. |2^, extend the analyses given 



in Ref. p4|, and will show that the present QMD approach gives a consistent description 



of the preequilibrium (p,xp') and (p,xn) angular distributions of ^''Al, ^^Ni and ^°Zr targets 
at 90 to 200 MeV energy range in the entire angular region without any fitting procedure. 
Based on the excellent agreement with the data, we then proceed to study some of the open 
problems left in the preequilibrium processes. In section ||, we will give a brief explanation 



of the QMD to show the essential feature of our model. In section |lTTt we compare our results 
with experimental data and predictions of the FKK theory to find similarities and differences 
between these two theories. We will then give a further discussion on the angular distribution 
of the quasifree scattering (QFS), and the energy dependence of the MPE process. 

II. BRIEF EXPLANATION OF THE QUANTUM MOLECULAR DYNAMICS 

A. Equation of motion 

We start from representing each nucleon (denoted by a subscript i) by a Gaussian wave 
packet in both the coordinate and momentum spaces in the following way: 

(r-R,)2 2L(p-Pi)=^" 



/i(r,p) = 8 ■ exp 



(1) 



AL 

where L is a parameter which represents the spacial spread of a wave paccket, Rj and Pj 
corresponding to the centers of a wave packet in the coordinate and momentum spaces, 
respectively. The total one-body phase-space distribution function is taken to be simply a 
sum of these single-particle wave packets. The equation of motion of Rj and Pj is given, on 
the basis of the time-dependent variational principle, by the Newtonian equation: 



R,: 



dH 



dH 



(2) 



and the stochastic N-N colhsion term as will be described below. We have adopted the 
Hamiltonian H to consist of the relativistic kinetic+mass energy and the Skyrme-type ef- 
fective N-N interaction p5| plus Coulomb and symmetry energy terms: 



i 



>' 



R-i — R.i 



erf I Ri -Rj I /V4L 



+ 7^7 E (1 - 2 I Q - 9 \)pij, 



(3) 



where "erf denotes the error function, and the Cj is 1 for proton and for neutron. The 
other symbols in eq.(|D are defined as: 

dp 

(4) 



= (27rL)-3/2exp [-(r - Ri)V2L 



and 



= E (47rL)-3/2exp [-(R, - R,)V4L 



(5) 



The symmetry energy coefficient Cg is taken to be 25 MeV. The four remaining parameters, 
the saturation density po, Skyrme parameters A, B and r are chosen to be 0.168 fm^^, —124 
MeV, 70.5 MeV and 4/3, respectively. These values give the binding energy /nucleon of 16 
MeV at the saturation density po and the compressibility of 237.7 MeV (soft equation of 
state, EOS) for nuclear matter limit. The only arbitrary parameter in QMD, i.e., the width 
parameter L, is fixed to be 2 fm^ to give stable ground state of taget nuclei in a wide mass 



range. These values are taken from our previous paper . 
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B. The collision term and the Pauli blocking 



The stochastic nucleon-nucleon coUision is taken into consideration as similar to the 
cascade model [23]: When the impact parameter of two nucleons is smaller than a value of 



(t/vt where a denotes the energy-dependent N-N cross section, an elastic or inelastic N-N 



collision takes place. We adopt a parametrization of N-N cross sections [23] which is similar 
to that of Cugnon [^ to take account of the in-medium effects which reduces the absolute 
magnitude and forward-peaking of the N-N cross sections. The angular distribution of the 
elastic scattering was selected by the Monte-Carlo sampling method. 

The Pauli blocking of the final phase-space is checked after each collision. The block- 
ing probability is calculated in the same way as the collision term in the Vlasov-Uehling- 
Uhlenbeck theory . 



The parameters in the N-N cross sections were fixed in Ref. [^ and are used in this paper. 
Together with the parameters of the one-body dynamics given in the previous subsection, 
the equation-of-motion of the QMD is uniquely determined. 

C. The ground state 

The ground state of the target nucleus is generated by packing Rj and Pj randomly 
based on the Woods-Saxon type distribution in the coordinate space and corresponding 
local Thomas-Fermi approximation in the momentum space, seeking a configuration to re- 
produce the binding energy calculated by the liquid-drop model within a certain (± O.SMeV) 
uncertainty. 

The average distribution of the Pj obtained for ^°Zr is shown in Fig. 1 as the solid 
histogram. Experimental nucleon momentum distribution in nuclei is parametrized by a 



superposition of 2-Gaussians [14 



as 
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where A^i is just a normalization constant, and the parameters po and go are related to the 
Fermi momentum p p via 



Po = \/'2/5pF, and 

go = \/3po- (7) 

This distribution is plotted as a broken curve in Fig. 1, where the parameter eo has been 
taken to be 0.07; about the mid-point of the range of this parameter recommended by 
Haneishi and Fujita. The nucleon momentum distribution in the QMD calculation has a 
similar shape to this 2-Gaussian distribution, while the commonly adopted uniform Fermi 
gas distribution is just a simple square-shaped function which vanishes above the Fermi 
momentum. The most significant difference among these distribution is the presence of 
the high-momentum component in the former two distributions which is not present in the 
uniform Fermi gas model. The presence of the high-momentum component is a common 
feature of finite- nucleon systems. As a matter of fact, the momentum distribution in QMD 
has a very similar shape to the one obtained by the Hartree-Fock theory as compared in 
Fig. 6(b) of Ref. [^. It is well known that the high momentum component enhances 
the backward angular distributions WM, and as will be shown later, we obtain the same 



conclusion from our QMD and 1-step Monte-Carlo simulations. However, the effect of the 
difference in the momentum distribution on the angular distribution of the primary particles 
emitted from the quasi-free scattering process was not very remarkable except at the very 
forward and backward angles, as will be discussed in the next section. It may worth noting 
here that the ground state in QMD as obtained in our work remains stable even with the 
high-momentum tail. 

D. Decomposition into step-wise contribution in multistep reactions 

For later discussion of the multistep reaction, it will be convenient to give a definition of 
step number in the QMD calculation which should reflect the number of collisions responsible 



7 



for emission of a nucleon. First we assign a step number of to each nucleon in the target 
nucleus. After a nucleon collides with incident nucleon, we set collision number 1 to both 
nucleons, inhibitting a collision between nucleons of collision number zero pair. Also, we 
prohibit successive collisions by the same partner. The rule of the change of the step number 
for each nucleon is that, if two nucleons i and j having step numbers Sj and Sj make a collision, 
the step numbers of both particles are modified to be Si + sj + 1. We then identify that a 
nucleon is emitted from the n-step process if an isolated nucleon emitted from the nucleus 
has a step number of n. 

As explained above, the first collision takes place only between the projectile and a 
nucleon in the target nucleus as expected intuitively. If one or two of these nucleons are 
emitted without experiencing further collisions, they contribute to the 1-step process. If, on 
the other hand, either of these nucleons makes a further collision in the nucleus, and one or 
two of these nucleons involved in the second collision are emitted without further collision, 
they are classified as the 2-step process. In the FKK theory, on the contrary, the 1-step direct 
cross section is calculated by means of the normal DWBA method averaged over many final 
Ip-lh states, that is caused as a result of having 1 coUision between the projectile and a 
nucleon in the target. The m-step FKK direct component is calculated by a folding integral 
of the (m-l)-step and 1-step cross sections, that results after a nucleon under interest has 
experienced the m-th collision in the system. The difinitions of the step number in QMD 
and FKK coincide up to the step number of 3. Beyond the 3-step process, however, those 
definitions become slightly inconsistent because QMD includes collisions between collided 
nucleons which are not present in FKK approach, although the probability of having such 
collisions in QMD is not very large. 



E. Calculation of the Cross Section 

In the calculation, many events having different impact parameter were generated. The 
impact parameter has been selected from a uniform distribution between and a maximum 
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J 2TTb ■ {M{E, n, b)) db (8) 



value which was taken to be shghtly bigger than the nuclear radius. The energy and direction 
of motion are stored event by event for every nucleon that becomes free (isolated from other 
nucleon), and finally the double-differential cross section was calculated as 

dE-dn 

where {M{E, Q, b)) denotes the average multiplicity of the particle under interest (neutron 
or proton) emitted in the unit energy-angular interval around E and Q for the impact 
paramter b event. 

Typically, 50000 events were generated to get a reasonable statistics in the step-wise 
double-differential cross section. In the calculation, the parameter has been fixed to the 



same values as in Ref. [23|, without any adjustment. 



III. RESULTS AND DISCUSSION 



A. Comparison with experimental data 



The calculated double-differential ^^Ni(p,xp') cross sections for incident energies at 120 
and 200 MeV, and the ^°Zr(p,xp') and (p ,xn) cross sections at 160 MeV are compared in Figs. 
2 and 3 with experimental data p8|-pO|] . The data have been shifted by the amount denoted 
in the parentheses. Agreement of the present calculation with the measured values is quite 
satisfactory from the very forward to backward angles, showing a basic ability and usefulness 
of our QMD approach to investigate the N-A reaction mechanisms in this energy regime. 
The problem of the underestimation at the backward angles in the semi-classical models 
PJT^ and the problems in the cascade model pT|-p3| at the very forward and backward 



angles are not present in the QMD approach. It must be also noticed that the QMD theory 
reproduces both the (p,xp') and (p,xn) cross sections simultaneously with a single set of 
parameters. This is a clear advantage of this approach over, e.g., the multistep direct FKK 
theory in which strength of the effective N-N cross section (the Vq parameter) must be 
adjusted depending on the projectile, ejectile, target and the incident energy Jl^. In this 
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way, it was verified that QMD gives an adjustment-free description of the pre-equihbrium 
(nucleon-in, nucleon-out) reactions at intermediate energy region in a unified manner. 



B. Step-wise contributions 



In order to have a better understanding of the reason why the QMD can reproduce 
the measured data so well, we compare in Fig. 4 the separate contributions to ^^Ni(p,xp') 
cross sections from the 1-, 2- and 3-step processes and the total of all steps calculated by the 
QMD theory with experimental data ||28|. Shown also are two arrows a and (5 corresponding 
to the angles expected from the 1-step quasifree scattering process without and with the 
acceleration effect by the mean field, i.e.. 



where i^out and E^-^ denote the energy of the outgoing and incoming particles in the laboratory 
frame, respectively, while V indicates the depth of the mean-field potential which has been 
taken to be -50 MeV. 

Fig. 4 indicates the followings: 

• the 1-step process is dominant at the forward angles, while at badkward angles the 2- 
and 3-steps are responsible to reproduce the measured cross sections. 

• the 1-step cross section does not have a peak neither at the angle a nor /?, instead it 
seems to have peaks at further forward angles for every secondary proton energy. As 
will be shown later, it is the Fermi motion of the target nucleon that is responsible for 
the shift of quasifree peak toward the forward angles. 

• the 1-step cross section does not fall off at the very forward angles for a high energy 
ejectile, i.e., at -Eout close to E-^^. This is a special feature of QMD theory, because 
the 1-step cross section calculated by the simple kinematical theory, as represented by 
the Kikuchi-Kawai formula , drops off sharply at the forward angles, which is the 
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reason why the cascade model often underpredicts the cross sections at this angular 
region. We will show later that the refraction of the projectile and the ejectile is 
responsible for not having the steep drop at the forward angles. 

• the 1-step cross section has non-negligible contribution beyond 90-degree. The mo- 
mentum distribution, especially the high-momentum component, is the reason of this 
spreading out the quasifree peak toward the backward angles. 

Therefore, three effects are found to be important to reproduce the measured (p,xp') cross 
section at the entire angular range; the refraction, the momentum distribution including the 
high-momentum component, and the multi-step contributions. Effects of the refraction and 
the momentum distribution will be discussed further in later subsections. 

C. Comparison with Feshbach-Kerman-Koonin model predictions 

The step-wise contributions to ^°Zr(p,xp') and (p,xn) reactions for incident energy at 
160 MeV, and ^''Al(p,xp') and (p,xn) reactions at 90 MeV calculated by the QMD are com- 
pared in Figs. 5, 6, 7 and 8 with those calculated by the multi-step direct FKK theory 
and the experimental data p9| , |30| , |35| , |36| . It is confirmed that both theories can reproduce 
the measured values rather satisfactorily. The similarity between the QMD and the FKK 
results, as well as their abilities to reproduce the data, are rather striking considering that 
these two theories are based on completely different concepts; the QMD being a superpo- 
sition of N-N scattering with mean-field effects, while the FKK is based on the DWBA 
scattering amplitudes. A noticeable difference, however, exists at the lowest ejectile energy 
of ^°Zr(p,xp') and (p,xn) reactions, where the FKK predictions are bigger than the measured 
data at forward angles {6 > 50-deg.), and are smaller at backward angles. The QMD results 
do not show such a problem. The main difference between the QMD and FKK calculations 
come from the difference in the 1-step cross sections; the 1-step FKK cross section has a 
prominent peak at around 30-degrees, and drops off steeply at backward angles, while the 
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1-step QMD cross section has much flatter shape. We will show later that the difference 
in the momentum distribution has little effect on the angular distribution shape from the 

1- step process. Therefore we conclude that the difference in the 1-step QMD and FKK cross 
sections come from the difference in the angular distribution of the elementary process; in 
the QMD calculation, the 1-step cross section is determined by the N-N cross section which 
is nearly isotropic in the CM of two colliding nucleons, while in the FKK theory it is de- 
termined by the DWBA. In spite of the difference in the 1-step cross sections, however, the 

2- and 3-step QMD and FKK angular distributions are very similar. This will be another 
conflrmation of the result obtained by Chadwick and Oblozinsky who have shown that 



the linear-momentum dependent state density obtained by the exact and statistical Gaus- 
sian solutions become identital at 2p-2h and 3p-3h states in spite of a difference in the Ip-lh 
state. 

D. Quasi- free scattering 

As shown in the previous sections, the 1-step quasifree scattering (QFS) cross sections 
calculated by the QMD theory has two prominent features; it does not fall off at the very for- 
ward angles unlike the kinematical calculations ||3^ , and it does not fall off at the backward 



angles as rapidly as one predicted by the FKK theory. Here, we investigate two items that 
may play important roles in the quasifree scattering process; the momentum distribution 
and the surface refraction effect. 

First, we have investigated the effect of the momentum distribution to QFS angular 
distribution. Figure 1 indicates that the momentum distribution in the QMD calculation 
differs noticeably from that of the uniform Fermi gas (UFG) model, which was adopted in 
Kikuchi-Kawai theory. Instead of the square-shaped distribution, the momentum distribu- 
tion in the QMD has a Gaussian-like shape with small portion above the Fermi momentum, 
which is in between the UFG and the 2-Gauss distribution suggested by Haneishi and Fujita 

0- 
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We have compared the QMD angular distribution from the 1-step (p,n) process of ^'^Zr for 
incident energy at 160 MeV in the topmost parts of Fig. 9 with a simple 1-step Monte-Carlo 
calculations with momentum distributions of UFG and 2-Gauss. The 1-step calculation was 
performed as follows: 

1. Select energy of a neutron in the target either from the UFG or from the 2-Gauss 
distributions, assuming a nucleon binding energy of 8MeV and Fermi energy (Ep) of 
40MeV. 

2. Make an isotropic scattering in the CM system of the projectile and the selected 
neutron in the target. 

3. The Pauli blocking effect is taken into consideration with a blocking probability given 

by 

Pblock = 1 - [1 - e{EF - E[)] ■ [1 - e{Ep - E'^)] (10) 

where E[ and £'2 denote the energies of scattered particles. 

4. If the collision is not blocked, the energy and angle of the scattered particle (which 
originally was in the target) in the laboratory frame is recorded. 

5. Repeat items 1 to 4 many times. 

6. The absolute magnitudes of these 1-step cross section was normalized to the corre- 
sponding 1-step QMD cross section. 

The upper two figures in Fig. 9 show that, in the main part of the angular distribu- 
tion, the difference between the UFG and the 2-Gauss results is not very noticeable. The 
main difference lies at the very forward and backward angles, where UFG result exhibits a 
steep drop, while the 2-Gauss result shows a slower decrease. This is definitely due to the 
high-momentum component in the 2-Gauss momentum distribution, because this difference 
disappears when we cut the high-momentum component in the 2-Gauss distribution . Also, 
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it is clear that the two distributions do not give the QFS peak at angles denoted by a nor 
/3, but give a peak at more forward angles. Therefore, the Fermi motion of target nucleons 
was found to shift the QFS peak to the forward angles. The peaks are, however, not very 
prominent in both cases; the Fermi motion tends to wash-out the QFS peak. The 1-step 
QMD cross section is in very good accord with both 1-step results at intermediate angles. 
As a matter of fact, the angle beyond which the QMD cross section vanishes lies in between 
the corresponding angles of the UFG and 2-Gauss results, because the momentum distri- 
bution in QMD lies in between these two distributions as shown in Fig. 1. However, the 
QMD results do not show decreasing angular shapes toward 0-degree. Therefore, the reason 
why the 1-step QMD results have large cross sections in the vicinity of 0-degree was not 
explained by the Fermi motion of the target nucleons. 

In the lower parts of Fig. 9, we have compared two kinds of 1-step QMD cross sections, 
one with the full calculation, and one which cuts the refraction effects. The 1-step QMD 
results without the refraction show a shape very similar to the one calculated with the 
uniform Fermi gas model, with a steep drop at the very forward angles, while the refraction 
effect totally washes out this steep decrease. Thus, it became clear now that it is the 
refraction effect by the mean-field which causes a non- decreasing 1-step cross sections in 
the QMD calculation at 0-degree region. This effect, together with contributions from the 
2-, 3- and higher steps, makes the total QMD cross sections to have a smoothly varying 
angular shape from the very forward to backward angles, which is in good accord with the 
measured data. Therefore, we conclude that the Fermi motion of target nucleons, mean field 
refraction, and the multi-step effects are essential in predicting the angular distributions of 
preequilibrium (N,N') cross section in this energy range. 

In the present calculation, the mean-field refraction effect washes out the decrease of the 
cross section at the very forward angle. This result, at first glance, may look completely 
opposite to the one obtained with the Geometry-Dependent Hybrid model (GDH) where the 
steep increase at the forward angle is washed out by the surface refraction at low incident 
energy (e.g.. Fig. 4 of Ref. At higher incident energy, however, the GDH predicts 
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decreasing angular shape toward 0-degree that is smeared out when the refraction effect is 
considered (20 MeV n of Fig. 6 in Ref. which is in good agreement with the present 
result. The reason why the GDH calculation has such a steep increase at the forward angle 
of low-incident energy is unclear to us: at least the 1-step process is concerned, the cross 
section must have a decrease as a result of kinematical restriction and Pauli-principle as 
Kikuchi-Kawai formula indicates. It may be just a result of finite angle-binning carried out 
by Blann et al. In the QMD calculation, both the Fermi motion and the surface refraction 
effect are taken into account by means of the equation-of-motion (i.e., Eq. in a unified 
manner, and any arbitrariness is not involved as is the case of GDH to introduce the surface 
refraction. 

It must be noticed that the refraction effect in our calculation might be overemphasized 
due to the fact that momentum dependence is not included in the effective N-N force that 
changes the mean-field potential from attractive at low energy to repulsive at energies higher 
than approximately 200 to 300 MeV region. However, the relativistic approach on the 
optical potential gives a wine-bottle-bottom shaped potential that remains attractive at 
the surface regime even at high energy region where the potential at the nuclear interior 
becomes significantly repulsive PE|J55[ ]. Moreover, according to Gadioli and Hodgson [|T|, 
inclusion of the momentum-dependent potential leads to two opposite consequences: (a) it 
reduces the importance of refractions of nucleons by reducing the potential at 200-300MeV 
region. This effect increases the probability of emission of particles. However, (b) the 
particles in the mean-field will have on the average lower kinetic energies, and (due to the 
increase of the N-N cross section) smaller mean free paths with the consequence that the 
nuclear transparency is decreased. These two effects tends to cancel each other, and the net 
effect of the momentum-dependent potential might be substantially reduced. 
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E. Multiple preequilibrium emission 



In Ref. [Q, we have shown that the QMD gives results consistent with the FKK theory 
for the energy spectra of first and second particles emitted from the preequilibrium process 
up to projectile energy of 200 MeV. Recently, there is a growing interest in nucleon- induced 
reactions up to 1 GeV region stimulated by the results obtained at LAMPF facility on the 
Pb(n,x) and Fe(p,x) reactions H^^Hl ^"^^ from a practical point of view P2| , ^ . Our QMD 
approach can be applied to above 1 GeV without any change, and there is no limitation 
in the number of particles emitted from the preequilibrium process. Therefore we have 
extended the analysis given in Ref. up to 1 GeV, and investigated the importance of 
MPE process as a function of projectile energy. 

We have shown in Fig. 10 the QMD results for the fractional contributions to the total 
number of inclusive proton emission contributed by the particles emitted as the 1st, 2nd, 
3rd, 4th and 5th particles. The upper figure shows the percentage when all emission energies 
are considered, while the lower figure includes only emissions of above 25 MeV (as defined in 
the preequilibrium regime in Ref. ||2^). The lower figure shows that sum of the contribution 
by the 3rd, 4th and 5th particles in the preequilibrium process occupies a fraction of about 
30 % at the incident energy of 500MeV, and almost 50 % at 1 GeV, showing a clear necessity 
of including the MPE process more than 2 particles. 



IV. SUMMARY 

We have shown that the Quantum Molecular Dynamics (QMD) model can reproduce 
the measured data for intermediate-energy nucleon-induced preequilibrium nucleon-emission 
process without any adjustment of the underlying parameters. Based on this success, we 
then have studied some of the open problems in the preequilbrium reactions; the angular 
distribution and the multiple preequilibrium particle emission (MPE). The QMD calculation 
has not shown the prominent quasifree scatterint peak, which is consistent with the measured 
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data. The reason of the overall agreement with the data was explained by the Fermi motion 
of target nucleons, the refraction of projectile and ejectile, and contribution from the multi- 
step processes. The MPE process beyond 2-particle emission was found to exceed 30 % at 
SOOMeV and reaches almost 50 % at 1 GeV, thus becoming the major reaction mechanism 
at this energy region. 
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FIGURES 

FIG. 1. Nucleon momentum distribution of ^°Zr. The solid histogram presents the resutls of 
QMD calculation. The smooth broken curve and the square distributions designates the 2-Gauss 
distribution parametrized by Haneishi and Fujita [15] and uniform Fermi gas distribution, respec- 
tively. 

FIG. 2. The 5^Ni(p,xp') cross sections at Ep = 120 MeV (left) and 200 MeV (right). The data 
have been multiplied by the amount denoted in the parentheses. 

FIG. 3. The ^'^Zr(p,xp') and (p,xn) cross sections at Ep = 160 MeV. The data have been 
multiplied by the amount denoted in the parentheses. 

FIG. 4. The 5SNi(p,xp') cross sections at Ep = 120 and 200 MeV. The total (thick solid line), 
1-step (dashed line), 2-step (broken line) and 3-step (thin solid line) QMD cross sections are 
compared with experimental data. The arrows a and /3 denote the position of QFS peaks as given 
by Eq. (9). 

FIG. 5. The ^°Zr(p,xp') cross sections at Ep = 160 MeV. The left row compares the predictions 
of total (thick solid line), 1-step (dashed line), 2-step (broken line) and 3-step (thin solid line) QMD 
cross sections with experimental data, while the right one those of FKK theory with experimental 
data. 

FIG. 6. The 9°Zr(p,xn) cross sections at Ep = 160 MeV. The left row compares the predictions 
of total (thick solid line), 1-step (dashed line), 2-step (broken line) and 3-step (thin solid line) QMD 
cross sections with experimental data, while the right one those of FKK theory with experimental 
data. 

FIG. 7. The 27Al(p,xp') cross sections at Ep = 90 MeV. The left row compares the predictions 
of total (thick solid line), 1-step (dashed line), 2-step (broken line) and 3-step (thin solid line) QMD 
cross sections with experimental data, while the right one those of FKK theory with experimental 
data. 
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FIG. 8. The ^^Al(p,xn) cross sections at Ep = 90 MeV. The left row compares the predictions 
of total (thick solid line), 1-step (dashed line), 2-step (broken line) and 3-step (thin solid line) QMD 
cross sections with experimental data, while the right one those of FKK theory with experimental 
data. 

FIG. 9. The 90Zr(p,xn) cross sections at Ep = 160 MeV and E^ = 120 (left figures) and En 
= 40 (right )MeV. The upper two figures compare 1-step QMD cross section (solid lines) with 
a simple 1-step Montc-Carlo calculations assuming the uniform Fermi gas (broken histograms) 
and 2-Gauss (long-broken histograms) momentum distributions of target nucleons. The lower two 
figures compare the full 1-step QMD cross section with a calculation ignoring the refraction effects. 

FIG. 10. Relative contributions of multi-particle emissions as a function of the incident energy 
calculated by the QMD. The upper figure includes all emission energies, whereas the lower figure 
considers only emission energies above 25 MeV (the preequilibrium regime). 
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